Theory for phonon-induced superconductivity in MgB 2 
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We analyze superonductivity in MgB2 observed below T c — 39 K resulting from electron-phonon 
coupling involving a mode at fuoi = 24 meV and most importantly the in-plane B-B E^g vibration 
at Tiui2 = 67 meV. The quasiparticles originating from tt- and cr-states couple strongly to the low- 
frequency mode and the i?2 9 -vibrations respectively. Using two-band Eliashberg theory, — 1.4 
and A CT = 0.7, we calculate the gap functions A* (a;, 0) (i = tt, a). Our results provide an explanation 
of recent tunneling experiments. We get Hfg/H^ fa 3.9. 
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It is important to understand the physics of the re- 
cently observed superconductivity in MgB 2 at T c «39K 
0. Recent experimental and theoretical results 
have shown that most likely strong electron-phonon in- 
teraction is the relevant mechanism for superconductiv- 
ity in this compound. However, many questions remain. 
In recent Raman ||, tunneling photoemission ||, 

and transport |hJ studies it was noted that the experi- 
mental data are compatible with a multi-component su- 
perconducting gap. This is consistent with the large 
observed anisotropy of H C2 |Tl| , |l2| ]. Inelastic neutron 
scattering(INS) experiments ]13| and point-contact spec- 
troscopy studies ]l4j reveal a peak at approximately 24 
meV in the phonon density of states in contrast to first 
principles calculations |15 17| . Recently, this peak was 
also indirectly confirmed by the analysis of the temper- 
ature dependence of the upper critical field H C2 (T) |jl8) . 
The isostructural AIB2 is not superconducting. 

In this paper we analyze superconductivity in MgB2 
using Eliashberg theory assuming that it arises due to 
electron-phonon coupling of the holes in the Boron planes 
(cr-states) to the E2 S -mode of B-B vibrations and due to 
the low-frequency mode interacting with the 7r-states. 

In Fig. 1 we illustrate the simple physical picture for 
superconductivity and the structure of MgB2. Similar 
to graphite, it has a quasi 2D-like structure, but with 
strong covalent bonds. According to previous studies 
j|,[||ll| , superconductivity occurs mainly in the B-planes 
involving two-dimensional cr- and and three-dimensional 
7r-states. In the first case the electron-phonon coupling 
can be estimated from g e - P h ^< <Jb\VV\<tb > • e, where 
the polarization vector e refers to the E2 9 -mode and 
ob refers essentially to the in-plane wave function of 
Boron bonds. Note, we obtain g e ~ P h ~ Vf ■ e (t denotes 
the corresponding hopping matrix element) and expect a 
resulting moderate electron-phonon coupling like in tran- 
sition metals. 

It is known from first principle calculations that MgB 2 
consists of four important phonon modes fl5| ]. Two 
modes involve vibrations of the Mg-B planes: the E\ u 
mode (40meV) corresponds to a sliding mode in xy- 



direction, and the Ai u mode (49meV) corresponds to 
a vibration in z-direction. On the other hand, the 
other two modes involve the B-B bonds: the E2 g mode 
(67meV) mentioned earlier corresponds to an in-plane 
breathing mode, and the B\ g mode (87meV) refers to an 
out-of-plane tilting mode of Boron atoms. First-principle 
calculations of the electron-phonon coupling o?F{lS) in 
Ref. jljl obtain that the dominating part of the pair- 
ing interaction stems from the Ei g mode. In the non- 
superconducting AIB2, this i?2g-mode lies at fiuj = 120 
meV. Recent INS and Raman experiments in MgB2 re- 
veal an additional peak in the phonon density of states 
at approximately 17-25meV 
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19[ | and, most impor- 
tantly, its softening in the superconducting state. It was 
noted [Q that this mode very likely is involving Mg vi- 
brations. Thus, one expects the coupling of this mode 
mainly to the p z - (tt) band due to its dispersion along the 
z-axis. Therefore, we conclude that superconductivity in 
MgB 2 is mainly induced by Boron tt- and cr-states cou- 
pled to low-frequency mode at 24meV and the i?2g-mode 
at 67meV respectively. Due to the different character of 
both states and their coupling to different phonon modes 
one expects the occurrence of superconducting gaps re- 
ferring mainly to the tt- and cr-bands. This is necessary 
as well for the analysis of the observed anisotropy in the 
upper critical field. 

For our analysis we use a two-band extension of 
Migdal-Eliashberg theory [ po|j2l| . We approximate the 
phonon spectral density F % (uj) by two Lorentzians cen- 
tered at fwJ\ ~ 24 meV for the 7r-band and Tw^-E^g) = 
67 meV for the cr-band. Hence, 
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describing the pronounced phonon peaks at the observed 
frequencies. Here, g p refers to the strength of the 
phonon mode i, and corresponds to its damping. For 
simplicity we neglect the q-dependence of the spectral 
functions, i.e. i*i(q) = 1. We employ Ti — 5 meV and 
T2 = 10 meV. Thus, the Eliashberg coupling constants 
(for a- and 7r-states) 
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Note, the single T c is provided by the non-zero value 
of A^cr and X aT! due to transitions between a- and 7r- 
states [^0|. However, their effect at low temperatures 
on the solution of the Eliashbcrg equations seems not 
to be of significant importance. Thus we neglect them 
for T=0. Regarding the value of X a we have chosen it 
0.7 in quantitative agreement with first principle calcula- 
tions . The strength of the coupling of the 7r-electrons 
to the low-frequency mode is unknown. However, us- 
ing an Einstein model for the phonon spectrum yielding 
a 2 F{uj) — g 2 /u>E, one expects A^ > X a [^2| (see also 
p3||). Taking into account also recent experiments indi- 
cating that the larger gap shows much smaller coherence 
peaks in the tunneling density of states |?J one concludes 
also that A^ > A CT . In the following, we approximate 
Att = 1.4. Then we use Eliashberg equations referring 
to the 7T- and a-states separately for T=0 (ui c is a cutoff 
frequency) : 
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where i refers to n- and a— electrons. Here Z 1 s (uj) de- 
notes the effective mass renormalization function of the 
a- and 7r-bands in the superconducting state and U % c 
refers to the screened Coulomb (pseudo-)potential calcu- 
lated within the static Thomas-Fermi approximation for 
3D metals. Solving these equations self-consistently for 
the given phonon spectrum (7r-states couple to u>i, while 
cr-states to w 2 -mode) we obtain A l (w,T) with structures 
due to the two pronounced modes in a 2 F l (uj) [Q. Re- 
garding the choice of the values for /i*, note that the 
a— band is more metallic than the 7r-band, which is re- 
flected by their different Fermi velocity and density of 



states at the Fermi level Jl|,|5|. This leads to n* a > fx*. 
Taking /i* = 0.1 and the appropriate values of TV^ep) 
and Wj, we obtain fi* = 0.18. These values of fi* and 
Xi give the observed magnitude of the superconducting 
gaps, their ratio and the anisotropy of the upper critical 
field H C 2- From this point of view, MgB2 can therefore be 
referred to as a conventional electron-phonon mediated 
superconductor with moderate coupling. 

In Fig.^ we show the results for the superconducting 
gaps for the a- and 7r-bands solving Eqs. ([|)-(^]). Due to 
the larger A and smaller /i* the gap value Aq = 7.5 meV 
for the 7r-band is almost two times larger than supercon- 
ducting gap Aq = 4.2 meV for the c-band. The values of 
the gaps and their ratio are in good agreement with ex- 
perimental data |^-|^Jlfj|] . The curves also reflect the un- 
derlying phonon spectrum a 2 F l {u). Re A l (w) changes 
its sign above uii due to the over-screened pairing po- 
tential and Im A l (w) starts to increase at uji because of 
phonon emission processes. 

In order to test whether the effective one-band model 
can be applied, we also have solved the one-band version 



of the Eliashberg equation for the average A 
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Att « 0.87 and average /2* = 0.5 (^i* + ^* ). The cor- 
responding phonon density of states includes both low- 
frequency and Ei g peaks with the coupling to the effec- 
tive Bloch states. The values of </; were chosen in order 
to have the same relative strength as in the case of the 
two-band model and also to satisfy the absolute value of 
A. Then, we find a superconducting gap Ao = 4 meV. 

In Fig. |^(a) we present our results for the tunnel- 
ing density of states for the effective gap. The obtained 
coherence peaks at 4 meV due to the effective supercon- 
ducting gap Ao agree with the experimentally observed 
position of the order of 4 meV 0. However, no fur- 
ther pronounced structure at approximately 7 meV as 
observed in experiment is obtained. This suggests that 
an effective one gap approximation is not a good one. 

In Fig.||(b), we show results for the two-band Eliash- 
berg theory. The tunneling density of states resulting 
from 7r- and a-bands is then given by 
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The results depicted in Fig.||(b) clearly show larger and 
smaller coherence peaks reflecting the corresponding su- 
perconducting gaps in the n- and c-bands respectively. 
The position of these peaks(AQ = 4.2 meV and AJ = 7.5 
meV) are in good agreement with experimental results 
(4.2 meV and 7.5 meV) indicating that these peaks may 
result from two superconducting gaps. The smaller co- 
herence peak of the larger gap is due to the smaller den- 
sity of states of the 7r-band at the Fermi level with respect 
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to the cr-states (see Eq. (0)). This is another indication 
that the larger gap has to be attributed to the 7r-states 
rather than to a er-band. 

In order to support further the validity of our re- 
sults, let us now turn to the discussion of the resulting 
anisotropy of the upper critical field H C 2. According to 
the anisotropic Ginzburg-Landau theory equations, the 
upper critical field along the c-axis is given by H£ 2 = ^p- 
and for the magnetic field applied in the a6-plane by 
~ TrfTbfT' Here, is the superconducting flux quan- 
tum, £ a ;, and £ c are the coherence lengths referring to the 
ab-plane and the c-axis in the clean limit, respectively. 
The coherence length is determined as £ w Thus, 
one gets the anisotropy ratio 
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Since only the 7r-band has a substantial dispersion along 
the fc z -axis, we can approximate vp by v F « 4 x 10 7 cm/s 
for the 7r-band pi. Similar arguments for the a— band 
lead to vf w vym 8.5 x 10 7 cm/s |j. Substituting our 

H ab 

values for the superconducting gaps we obtain j^- m 3.9, 

c 2 

which is in good agreement with the experimental obser- 
vation 11 The recent suggestion [§6],^7| that the 
larger gap belongs to the a— band and the smaller gap 
occurs in the 7r-band reduces significantly the anisotropy 
of the upper critical field and also leads to the intensity 
ratios of the coherence peaks in the tunneling density 
of states being inconsistent with experiment J7J. How- 
ever, we get H® 2 b ~ 5.6Tesla taking into account the 
renormalization of the Fermi velocities due to electron- 
phonon interaction, while some experiments give 15Tesla 
or even large values. Note, further correcting the LDA 
results for vp by a factor 1.5 for the it- and a— states we 
get w lA.QTesla. Many-body effects like electron- 
electron correlations may cause this correction. 

In FigJ| we show the calculated effective mass 
renormalization functions Z l s (uS) for the a- and 7T- 
quasiparticles reflecting many-body effects due to 
electron-phonon interaction and which describe the 
strength of the interaction between the corresponding 
quasiparticles. Their structures again reflect the under- 
lying phonon spectrum F l (ui) at 24 meV for the 7r-states 
and at 67 meV for the cr-electrons. Remarkably, the mass 
renormalization for 7r-electrons is larger due the differ- 
ence in the corresponding A^. 

Let us now estimate the high superconducting tran- 
sition temperature and the corresponding isotope effect. 
As we already mentioned the single T c is expected due to 
non-zero values of A„. Moreover, as was shown earlier 
p0| , even in a two-band case one still can use effective 
parameters for the determination of T c . Thus, we use 
the McMillan formula (w 2 ~ 1/VMb) PI 
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1.04(1 + A) 
A-/2*(1 + 0.62A) 
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with uip> uj2, and an effective A = 0.87 and Jx* = 0.14 
resulting in T c as 32K. We get for isotope coefficient 
a = — dh\T c / dhiMs the value 0.43. Assuming that the 
isotope effect will be mainly determined by the coupling 
between cr-states and the E2 g -mode we should use in Eq. 
(9) for A and fi* the values referring to the cr-states. Thus, 
we find the absolute value a «0.23 which is in agreement 
with experiment . We also would like to point out that 
despite of the larger gap in the 7r-band, superconductiv- 
ity (n s ) and T c are mainly due to the cr— states. This 
can be seen from the estimation of the superfluid density 
n\ ~ Ni(ep)u>i/(xmit volume) which is much larger for 
the cr-electrons due to their larger density of states at 
the Fermi level and the large u>2- 

In summary, we can explain important facts about 
superconductivity in MgB2 as resulting from electron- 
phonon coupling of the 7r-states and the cr-states with 
the corresponding phonon modes. The larger gap results 
from the 7r-band coupled strongly to the low-frequency 
mode, while the smaller gap is due to the coupling of cr- 
bands with the E2 g -mode. The structure seen by recent 
tunneling experiments is due to the larger superconduct- 
ing gap of the 7r-band. We also find H£/HZ 2 w 3.9 
in good agreement with experiment. We safely conclude 
that there is no need for invoking other pairing mecha- 
nisms besides electron-phonon interaction. Direct calcu- 
lations of 3e/-pfc an d \ are necessary for a further sup- 
port of our simple model. In particular, the origin of 
the 24 meV peak observed by INS and Raman scattering 
must be clarified. Using our picture for superconductiv- 
ity in MgB2, one may estimate the absence of supercon- 
ductivity in AIB2 due to the decrease of X a (A ~ \g\ 2 /0J2) 
and the shift of cr-states below ep Q . 
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FIG. 1. Simplified physical picture for superconductivity 
in MgB2: a) Top view of the elementary unit cell of MgB2 
consisting of hexagonal layers of Mg and graphite-like Boron 
layer. The Boron a-states (p x and p y ) are involved in super- 
conductivity and couple to the Ei g mode of B-B vibrations 
shown by arrows; b) the Boron 7r-states (p 2 ) are coupled to 
the out-of-plane vibrations shown schematically by arrows. 



FIG. 2. Results for the gap function A^u) at T = OK 
using two-band Eliashberg-theory: (a) cr-band with Ao- = 0.7 
and f£ = 0.18 yielding AJ = A CT (w = 0) = 4.2 meV. (b) 
7r-band with A^ = 1.4 and fj,* = 0.1 yielding AJ = 7.5 meV. 



FIG. 3. Calculated tunneling density of states: (a) using 
an effective one-band Eliashberg theory with elastic scattering 
strength 7 = 0.3meV; (b) using two-band Eliashberg theory 
and the same 7 as in (a). In the inset, experimental data is 
displayed M. 



FIG. 4. Results for the mass renormalization function 
Zl(u>) at T = OK using two-band Eliashberg-theory for (a) 
the (T-band and (b) the 7r-band. 
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